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Abstract 

We study the renormalization group evolution and the infra-red stable fixed 
points of the Yukawa couplings of the non-minimal supersymmetric standard model 
(NMSSM) with R-parity violation. Retaining only the R-parity violating couplings 
of higher generations, we analytically study the solutions of the renormalization 
group equations of all the couplings of the model. We find that there are no si- 
multaneous non-trivial infra-red fixed points for all the couplings of the model, and 
that the infra-red fixed point structure of the model is similar to the MSSM with 
R-parity violation. In particular, we show that only the baryon number violating 
coupling A233, together with top- and bottom-quark Yukawa couplings, approaches 
a non-trivial infra-red stable fixed point. However, this fixed point solution predicts 
a top-quark Yukawa coupling which is incompatible with the top quark mass for 
any value of tan j3. 
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Recently considerable attention has been devoted to the study of infra-red (IR) sta- 
ble fixed points of the Yukawa couplings of the minimal supersymmetric standard 
model 0. In the minimal supersymmetric standard model(MSSM) 0, gauge invariance 
and supersymmetry allow dimension four couplings which violate baryon (B) and 
lepton number (L). These couplings are usually forbidden by invoking a discrete sym- 
metry 1^ called R-parity {Rp)- However, the assmption of Rp conservation at the level 
of MSSM appears to be ad hoc, since it is not required for the internal consistency of 
the model. Much attention has, therefore, been paid to the study of the renormalization 
group (RG) evolution of these B- and L-violating (and Rp violating) Yukawa couplings. 
This includes the study of the quasi-fixed point behaviour ^ , as well as the true infra-red 
fixed points of the different Yukawa couplings and the analysis of their stability 0. This 
has led to certain insights and constraints on the fixed point behavior of some of the 
Rp violating Yukawa couplings, involving higher generation indices. For the fixed point 
or quasi-fixed point scenarios to be successful, it is necessary that these fixed points be 
stable [^]. It has been shown that in the minimal supersymmetric standard model with 
R-parity violation, only the S- violating coupling A233, together with top- and bottom- 
quark Yukawa couplings, approaches a non-trivial infra-red stable fixed point, whereas all 
other non-trivial fixed point solutions are either unphysical or unstable in the infra-red 
region 0. 

It is well known that the minimal supersymmetric standard model suffers from the 
so-called fi problem associated with the bilinear term connecting the two-Higgs doublet 
superfields in the superpotential. A simple solution to this problem is to postulate the 
existence of a Higgs singlet superfield, and to couple it to the two Higgs doublets in the 
superpotential via a dimensionless trilinear coupling. When the Higgs singlet obtains a 
vacuum expectation value, a bilinear term involving the two Higgs doublets is naturally 
generated[P]. However, this leads to additional trilinear superpotential couplings in the 
model, the so called non-minimal supersymmetric standard model(NMSSM). Because of 
the additional trilinear couplings, it is important to study the infra-red fixed points of 
the NMSSM, and analyze the effect of these additional trilinear Yukawa couplings on the 
infra-red behaviour the other Yukawa and B- and L-violating couplings, and contrast it 
with the situation that obtains in the MSSM. 

To this end we shall consider renormalization group evolution and infra-red fixed points 
of the Yukawa couplings of the non-minimal supersymmetric standard model (NMSSM) 
with baryon and lepton number (and Rp) violation. The superpotential of the NMSSM 
is given by 

W = {hu)abQiURH, + {hn)abQiD'^H^ + {hE)abLlE'^H, + XNH^H^ - (1) 

Gauge invariance, supersymmetry and renormalizability allow the following superpotential 
B and L violating terms in the NMSSM: 

Wb = ^x'LdId'^u% (2) 
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where the notation is standard, except that there is a new L-violating term with the 
dimensionless Yukawa couphng Aa in (^. This term can be rotated away into the R- 
parity conserving term XNH1H2 via an SU{A) rotation between the superfields Hi and 
Li. However, this rotation must be performed at some energy scale, and the term is 
regenerated through the renormahzation group equations. The Yukawa couphngs Xabc 
and X'^fj^ are antisymmetric in their first two indices due to SU{2)l and SU{?>)c group 
structure, respectively. Phenomenological studies of supersymmetric models of this type 
have placed constraints [|10] on the various couplings Xahc, A^^^ and A"^^ in the MSSM, but 



there is still considerable room left. We note that the simultaneous presence of the terms 
in Eq. @j and Eq. (|) is essentially ruled out by the stringent constraints [|ll| implied by 
the lack of observation of nucleon decay. 

In addition to the dominant third generation Yukawa couplings ht = (/i[/)33, hf, = 
(^0)33 and hj. = (/i£)33, and the trilinear couplings A and k in the superpotential (|1]), 
there are 9 independent Rp violating couplings A^'^^ in Eq. (H), and 39 independent Xabc, 
X'^Ij^ and A^ in Eq. (^). Thus, one would have to solve 14 coupled non-linear evolution 
equations in the B-violating case, and 44 in the L-violating case, in order to study the 
evolution of the Yukawa couplings in the nonminimal supersymmetric standard model 
with Rp violation. In order to render the Yukawa coupling evolution tractable, we need 
to make certain plausible simplifications [0]. As is usually done in the analysis of 
the MSSM with R-parity violation, we shall retain only the highest generation couplings 
A3, A233, A333 and A233, and neglect the rest. We note that the Rp violating couplings 
to higher generations evolve more strongly because of larger Higgs couplings in their 
evolution equations, and hence could take larger values than the corresponding couplings 
to the lighter generations. Furthermore, the experimental upper limits are stronger for 
the Rp violating Yukawa couplings corresponding to the lighter generations. 

We shall first consider the RG evolution of Yukawa couplings arising from superpo- 
tential terms in ([^) and (0), which involve baryon number violation. We have derived the 
one-loop renormahzation group equations for all the trilinear couplings of the model. For 
the couplings ht, hb, hr, A, k and A233 (all other Yukawa couplings set to zero) these can 
be written as: 

2 dht _ , 2 , . 2 , \2 , o\"2 16 2 o„2 13 2 



16>r=5^ = h {hi + 6hl + ft? + A' + 2A« - fgl - 3gl - ^g? 

2 dhr _ , /o!,2 , ^7,2 , x2 o„2 9 n 



levr^^^ = A (3/^? + 3hl + hl + 4X' + 2e - 3gl - ^1 
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We note that evolution equations for the gauge couphngs {gi,i = 1,2,3) in NMSSM are 
identical to those in the MSSM, since the additional Yukawa couplings do not play a role 
at this order. With the definitions 
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(5) 



ordering the ratios as Ri = {R", Rr, Rt, Rt, Rx, Rk), and retaining only the SU (3)c gauge 
coupling, we can rewrite the renormalization group equations @ in the form (t = — /n/i^) 



dRi 

~dr 



asRi 



(ri + &3) - X! ^ij^j 



(6) 



where = J2r 2 Cr, Cr is the QCD Casimir for the various fields {Cq = Cjj = Cjy = 4/3), 
the sum is over the representation of the three fields associated with the trilinear coupling 
that enters Ri, and is a matrix whose value is fully specified by the wavefunction 
anomalous dimensions. A fixed point, denoted as R*, is then reached when the right 
hand side of Eq. (||) is for all i. We, then, see that there are two fixed point values for 
each coupling: R* = 0, or the non-trivial fixed point solution 



R: = Y.is-%ir,+b,). 



(7) 



Since we shall consider the fixed points of the couplings ht, A, k and A233 only, we shall 
ignore the evolution equation for h^-. However, the coupling h^- does enter the evolution 
of hb, but it can be related to hb at the weak scale (which we take to be the top-quark 
mass), since 

\/2mT-{nit) mrirrir) 7]b 



K{mt) 



-hb{mt) = 0.6hb{mt), 



rjrV cos (3 mb{mb) rjr 

where rjb gives the QCD or QED running [|12| of the b-quark mass mb{^) between fi = rrib 
and jj, = nit (similarly for rfr), and tan [3 =< H2 > / < >, with v = {^/2Gf) 
GeV. The anomalous dimension matrix S can, then, be written as 



'1/2 



246 



S 



( 6 
2 

2 1 
3 + r/ 



\ 



2 2 \ 

6 + r/ 1 1 

6 1 

3 4 2 

6 6 / 



(9) 



with ordering {R" , Rb, Rt, Rx, Rk) , where t] = hl{mt) / hl{mt) ~ 0.4 is the factor coming 
from Eq. (H). We, therefore, obtain the following fixed point solution for the ratios: 



/ p//* p* D* D* D*^ / 



25 



5(5 + 7]) 115 + 2477 95 + 227] 



12' 4(10 + ?])' 4(10 + r/)' 4(10 + r/) ' 4(10 + r/) 



(10) 
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Since R\ is negative, this is not an acceptable fixed point solution. Thus, there is no 
infra-red fixed point solution for NMSSM with baryon number violation where all the 
trilinear superpotential couplings attain non-zero fixed point values. 

We next try to find a fixed point solution with Rl = 0, and R"*, R^, Rl and Rl being 
given by their non-zero solutions, which is relevant for the low values of the parameter 
tan j3. We need to consider the appropriate 4x4 sub-matrix of the matrix S in Eq. (^) 
to obtain the fixed point solutions for R"*, R^, R^ and Rl in this case. This sub- matrix is 
given by 

^ 6 2 
2 6 10 
3 4 2 

V 6 6 y 



S 



'ID 



with the fixed point solution 



Rl = 0, 



. 95 10 
138' 23^ 



38 53, 
23' 46' 



(12) 



This is not a theoretically acceptable solution, as all the fixed point values are not positive. 

We now try to find an infra-red fixed point solution with the trlinear coupling R^ = 0. 
Proceeding in the same manner as in the case of Rl 

. 385 + 76r/ 



RX =0, 



{R" 



p* p* D*] 



0, we get the fixed point solution 
44 10 + 2r] 



12(85 + 167]) ' 3(85 + 16r/) ' 85 + 16r/' 



1. 
'2' 



(13) 



Since Rl < this fixed point solution must be rejected. We can also try to find a fixed 
point solution with Rl = with the result 



Rl = 0, 



{R" 



p* p* p*\ 



, 455 + 83r/ 37 51(5 + r/) _ 435 + 96?] 

'30(25 + 4?])' 6(25 + 47])' 30(25 + 477) '~ 30(25 + 477) ' 



(14) 

which is not an acceptable solution either. Finally, we try a fixed point solution with 
R"* = 0, and all other Yukawa couplings attaining non-zero fixed point values. This is 
the case of NMSSM with R-parity conservation. In this case we find the solution 



R" 



0, 



(Rl 



P* p* p*"! 



5{5 + ri) 105 + 23r7 90 + 19r7 



'3(10 + r]) ' 3(10 + 7]) ' 3(10 + r]) ' 6(10 + r]) ' 



(15) 



which must also be rejected. We have, thus, shown that there is no infra-red fixed fixed 
solution with one of the trilinear couplings being zero, and all others attaining a nonzero 
fixed point value. 

Having failed to find an acceptable fixed point solution with one of the couplings being 
zero, we now try to find a solution with two of the trilinear coupling attaining a zero fixed 
point value. In this case we find only the following acceptable fixed point solution: 



p* p* 



0, 



{R" 



p* p*^ 



, 385 + 76r7 10 10 + 2r/) , 

'3(170 + 32r/)' 85 + I677' 85 + 167?' 



(16) 
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Interestingly, these fixed point values for the couplings R"*,Rl and are same as in 
the minimal supersymmetric standard model with third generation baryon number viola- 
tion 0. 

We can also try to find fixed point solutions where three of the trilinear couplings 
attain zero fixed point values, whereas the remaining two attain non-trivial fixed point 
values. In this case we find the following theoretically acceptable fixed points: 

19 1 

Rl = Rl = Rl = 0, iR"\ r;) = g), (17) 

35 35 _|_ 7^ 

Rl = Rl = R"* = 0, (Rl Rl) = (3(35 ^g^), 3(35 + 6r^))' (18) 

2 7R _|_ 1 5r) 

fl: = fi: = fi,- = o. (is;. = (_, g^^). (19) 



All the infra-red fixed point solutions in {\17\} - (^) have one thing in common, namely 
the trilinear couplings A and k approach the fixed point value zero in the infra-red region. 
Furthermore, the other trilinear couplings approach the same infra-red fixed point values 
as in MSSM §. 

Having obtained more than one theoretically acceptable infra-red fixed points in the 
NMSSM with baryon number violation, it is important to determine which, if any, is likely 
to be realized in nature. To this end we must examine the stability of each of the fixed 
point solutions ([TBI) - (0). 

We first consider the stability of the fixed point solution (|16|). Since in this case the 
fixed points of the couplings R\ = Rl = 0, we have to obtain the behaviour of these 
couplings around the origin. This behaviour is determined by the eigenvalues S 



1 



5 

Y,S,,R*-{r, + h 
i=3 



1,2, (20) 



where ri = r2 = 0, and the matrix S is the matrix appearing in the corresponding RG 
equation (H) but now with the ordering of the ratios as Ri = {Rx, Rk, R", Rb, Rt), and the 
fixed points R*, i = 3,4,5 corresponding to the fixed point values of R"*, Rl, and R^ in 
Eq. (|16|). Inserting these values in Eq. (pO]), we find 

\X] -\ + 11 f21) 

thereby indicating that the fixed point is attractive in the infra-red direction. The be- 
haviour of the couplings R", Rb and Rt around their respective fixed points is governed 
by the sign of the eigenvalues of the matrix A whose entries are {i not summed over) p 

= 'T~R*i^ij^ (22) 

where R* = {R"*,Rl,Rl) is the fixed point solution in (0), and Sij is the upper left 
corner 3x3 submatrix of the matrix (0). For stability, we require all the eigenvalues of 
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the matrix (|2^) to have negative real parts (note that the QCD /3-function 63 = —3 is 



negative). The eigenvalues of the matrix ( p2D are calculated to be 

[A,].=3,4,5 = [-1-6,-0.2,-0.2] (23) 

which shows that the fixed point (p!6D is an infra-red stable fixed point. We note that 
the eigenvalue A3 is larger in magnitude as compared to the other eigenvalues in (^31) , 
indicating that the non-trivial fixed point for A233 is more attractive, and hence more 
relevant. 

Next, we consider the stability of the fixed point solution ([171). Since in this case the 
fixed point of the couplings = Rl = Rl = 0, we have to obtain the behaviour of these 
coupling around the origin. This behaviour is determined by the eigenvalue [pf 



^■4 



=4 



1,2,3, (24) 



where ri = r2 = 0, = 2{Cq + Cjj) = 16/3, and the matrix S is the matrix appearing in 
the corresponding RG equation (D with ordering of the ratios as Ri = {Rx, Rk, Rb, R", Rt), 
and the fixed points R*, i = 4, 5 corresponding to the fixed point values of R"* and R^ in 
Eq. (p!?!). Inserting these values in Eq. (0), we find 



[A.U,2,3 = [-^, -1, (25) 

showing thereby that the fixed point ( P!7D is unstable in the infra-red region. Similarly, it 
can be shown that the fixed point solutions (|18|) and ([19|) are unstable fixed points. 

One may also consider the case where the couplings X2^^,hh,X and k attain trivial 
fixed point values, whereas ht attains a non-trivial fixed point value. In this case we have 
RI = R^ = 7/18, which is the same as the Pendleton- Ross |l[] top-quark fixed point of 
the MSSM. We must, of course, study the stability of this solution in the present context. 
To do so, we must consider the eigenvalues 

A. = Us^5Rt - (r, + 63)), i = 1,2,3,4, (26) 
03 

where Si^ are read off from the matrix (P), but now with the couplings reorderd according 
to Ri = {R", Rb, Rx, Rk, Rt), thereby yielding 

_ 38 35 25 

[Aijj=l,2,3,4 - ~ ~ ^' ^ ' 

Since the sign of each of Ai and A2 is positive, this solution is also unstable in the infra-red 
region. Nevertheless, from our discussion of infra-red fixed point solution (0), it is clear 
that the Pendelton-Ross fixed point would be stable in the NMSSM in case hi, and A233 
are small, though negligible at the GUT scale (with X = k = 0). In this case, these would, 
of course, evolve away from zero at the weak scale, though realistically they would still 
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be small (but not zero) at the weak scale. Thus, the only true infra-red stable fixed point 
solution is the baryon number, and violating solution (|1^). We note that this fixed 
point solution is identical to the corrresponding fixed point solution in the MSSM with 
baryon number violation This is one of the main conclusions of this paper. We note 
that the value of in (|1^) is lower than the corresponding value of 7/18 in MSSM and 
NMSSM with Rp conservation. 

It is appropriate to examine the implications of the value of ht{mt) predicted by our 
fixed point analysis for the top-quark mass. From (|T6|), and a^{mt) — 0.1, the fixed point 
value for the top- Yukawa coupling is predicted to be ht{mt) ~ 0.4. This translates into 
a top-quark (pole) mass of about rrit ~ 70 sin /3 GeV, which is incompatible with the 
measured value of top mass, ~ 174 GeV, for any value of tan/?. It follows that 
the true fixed point obtained here provides only a qualitative understanding of the top 
quark mass in NMSSM with Rp violation. 

We now turn to the study of the renormalization group evolution for the lepton num- 
ber, and Rp, violating couplings in the superpotential (0). Here we shall consider the 
dimensionless couplings A3, A233 and A333 only. Furthermore, we shall restrict our atten- 
tion to one kind of lepton number violation at a time. Thus, we shall consider three 
different cases, i.e., we shall take A333 ^ A3, A 



or A233 ^ '^3 5A333, or A3 ^ A2335A333, 



respectively. In the case when A333 is the dominant of the lepton number violating cou- 
plings, we define R = X^ss/ds^ reorder the couplings as Ri = {R , Rf,, Rt, Rx, Rk), so 
that the RGEs for this case can be written as 

dR, 



dt 



(28) 



where S is the anomalous dimension matrix 



S 



6 6 + r] 1 

6 Q + T] 1 1 

1 1 6 10 

3 + r] 3 4 2 

6 6 



(29) 



1] = h'^{mt) / hl{mt) — 0.4 is the factor coming from Eq. (^, and other quantities are 
defined in a manner similar to the case of baryon number violation. The non-trivial fixed 
points are given by Eq.(|^), and are calculated to be 

, 420 + 92r] 315 105 + 29r/ 



/ p * p* p* p* p* 



0,-2). (30) 



'315 + 11477 315 -M14r/' 315 + 11477' 

Since all the fixed point values are not positive, this is an unphysical solution. Thus, we 
conclude that a simultaneous non-trivial fixed point for all the couplings A333, hf,, ht, A and 
k does not exist. We next try to find a fixed point solution with Rl = 0, and R'*, R^, i?^ 
and Rl being given by their non-zero solutions, relevant for low tan /3 values. In this case 
we find the fixed point solution 

R, = 0, {R ^Rt,Rx,Rk) = {^,Y8'-Y^,Y^), (31) 
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which contains a non-physical negative value, and hence is not an acceptable fixed point 
solution. We next try to find a solution with two of the trilinear couplings attaining a zero 
fixed point value. We found that there are no such physically acceptable fixed points. We 
then tried to find fixed points with three of the trilinear couplings attaining a zero fixed 
point value. There is only one such acceptable fixed point: Rl = = Rl = 0, R* = 
Rl = ^. However, stability analysis shows that this fixed point is either a saddle point, 
or an ultra-violet fixed point. Furthermore, there are no fixed points with four of the 
couplings attaining a zero fixed point value. We conclude that there are no non-trivial 
stable fixed points in the infra-red region for the lepton number violating coupling A233. 

If on the other hand the coupling A233 is the dominant of the lepton number violating 
couplings, we define R = Xlss/ di^ reorder the Yukawa couplings as {R, Rb, Rt, R\, Rk), 
so that the relevant RGBs can be written as 



dRi 



azRi[{ri + 63) - ^SijR 



'3\i 



(32) 



with the anomalous dimension matrix 
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(33) 



with the by now usual definitions of the various quantities in 
fixed point is then calculated to be 



and 



The non-trivial 



{R* 



TD* T3* TD* TD*^ 



90 + 109// 25 5(5 + r/) _105 + 23r/ 30 + 43// 
12(10 + 7])' 3(10 + r/)' 3(10 + 7])' ~ 3(10 + r/) ' 6(10 + 7/)' 



(34) 

which must be rejected as unphysical. We next try to find a fixed point solution with 
i?^ = 0, and -R'*, i?^ , R\ and R^. being given by their non-zero solutions. In this case we 
find the fixed point solution 



D* n ( T?* D* D* D*^ 



3 20 19 29, 
'4' 27'~"9'' 18' 



(35) 



which, again, is physically unacceptable. Proceeding further as in the previous cases, we 
have not found any physically acceptable stable infra-red fixed points in this case. Thus, 
there are no fixed point solutions for the lepton number violating coupling A233. 

Finally, if the lepton number violating coupling A3 is the dominant coupling, we define 
Rz = ^l/dh reorder the couplings as {R3, Rb, Rt, R\, Rk), so that the RGBs can be 
written in the form 

dR- 

asRi[{ri + 63) - X! ^ijRj]^ (36) 
j 



dt 
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with the anomalous dimension matrix 
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(37) 



which leads to the fixed point values for the couplings 

{ D* D* D* D* D* \ ( 



100 20 43 29 
"27' '27' 27' 18^' 



(38) 



This fixed point is physically unacceptable. Thus, in this case also there is no fixed point 
with all couplings approaching non-trivial fixed point values. 

We can next try a fixed point with = 0, which is relevant for low values of tan /3. 
In this case we find the fixed point 



20 
27' 



19 



29. 
18^ 



(39) 



which is again unacceptable. We now try to obtain fixed point solution with two of the 
couplings R^ = Rl = 0, and other couplings having non-trivial fixed point values. We 
find 

20 19 29 

Rl = Rl = 0, (i?*, Rl Rl) = (-, --,-), (40) 



V i' A) kj V27' 9 ' 1^ 
which must also be rejected. We have also tried to obtain a fixed point solution with 

RX 



D* 



0, with other couplings having non-zero fixed point values. We find 

210 + 55?7 55 97 + 22r7 , 



D* D* n ( E?* E?* E?*^ 



3(61 + llr^)' 3(61 + 11?])' 3(61 + Ur])' 



(41) 



which, again, is unacceptable. We have also checked that the trivial fixed point for the 
couplings A3, hb, A, k and the Peddelton-Ross type fixed point for the top-quark Yukawa 
couphng is unstable in the infra-red region. We, therefore, conclude that in this case of 
lepton number violation also there are no acceptable infra-red fixed points. 

Thus, there are no non-trivial, infra-red stable fixed points for any of the lepton number 
violating couplings in the non-minimal supersymmetric standard model. 

To conclude, we have analyzed the one-loop renormalization group equations for the 
evolution of Yukawa couplings in the non-minimal supersymmetric standard model with 
Rp violating couplings to the heaviest generation, taking into account B and L violating 
couplings one at a time. The analysis of model yields the surprising and important 
result that only the simultaneous non-trivial fixed point for the baryon number violating 
coupling A233 and the top-quark and b-quark Yukawa couplings ht and hb, and the trivial 
fixed point for A and k, is stable in the infra-red region. However, the fixed point value 
for the top-quark coupling here is lower than its corresponding value in the MSSM, and 
NMSSM, with R-parity conservation, and is incompatible with the measured value of the 
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top-quark mass. Thus, it appears that the baryon number, and Rp, violating couphng has 
the effect of reducing the infra-red fixed point value of the top-quark Yukawa coupling to 
the same extent in MSSM and NMSSM. The Rp conserving solution with A233 attaining its 
trivial fixed point, with ht and hh attaining non-trivial fixed points, is infra-red unstable, 
as is the case for trivial fixed points for A233 and hb, with a non-trivial fixed point for 
ht- Our analysis shows that the usual Pendleton- Ross type of infra-red fixed point is 
unstable in the presence of Rp violation, though for small, but negligible, values of hb and 
A233 it could be stable. We have also found that there are no non-trivial infra-red stable 
fixed points for the lepton-number, and R-parity, violating couplings in the NMSSM. Our 
results are the first in placing strong theoretical constraints on the nature of Rp violating 
couplings in the NMSSM from fixed-point and stability considerations: the fixed points 
that are unstable, or the fixed point that is a saddle point, cannot be realized in the 
infra-red region. The fixed points obtained in this work are the true fixed points, and 
serve a lower bound on the relevant Rp violating Yukawa couplings. Their structure is 
essentially the same as in MSSM with R-parity violation. In particular, from our analysis 
of the simultaneous (stable) fixed point for the baryon number violating coupling A233 and 
the top and bottom Yukawa couplings, we infer a lower bound on A233 ~ 0.98. 
Acknowledgements: The work of (PNP) is supported by the University Grants Commission 
research award No. 30-63/98/SA-III. 
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